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Abstract. This contribution presents a Bayesian approach to the issue of linking of the 
results from key comparison measurements. A mathematical treatment based on Bayesian 
statistics for the analysis of the results from two comparisons with some joint participants is 
described. This robust statistical analysis provides expressions and standard uncertainties for 
the key comparison reference value (KCRV) and the degree of equivalence (DOE) as well as a 
conformity check without any assumption on a priority of one of the comparisons. In addition 
to the derivation of the mathematical formulae to be used for this type of ‘distributed linking’, 
we also present one synthetic and one real linking example and discuss possible applications 
of this new linking procedure. 
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1. Introduction 

The CIPM Mutual Recognition Arrangement (MRA) [I], which was put into force in 
October 1999 and meanwhile was signed by 95 National Metrology Institutes (NMIs) and 
4 international organisations worldwide describes the foundations and requirements for the 
mutual recognition of national measurement standards and for calibration and measurement 
certificates issued by NMIs. In addition the CIPM MRA covers a further 150 institutes 
designated by the signatory bodies. The MRA, chapter 3 states, that “... the technical 
basis of this arrangement is the set of results obtained in the course of time through key 
comparisons carried out by the Consultative Committees of the CIPM, the BIPM, and the 
regional metrology organisations (RMOs), and published by the BIPM and maintained in the 
key comparison database. Key comparisons carried out by Consultative Committees or the 
BIPM are referred to as CIPM key comparisons; key comparisons carried out by regional 
metrology organisations are referred to as RMO key comparisons; RMO key comparisons 
must be linked to the corresponding CIPM key comparisons by means of joint participants. 
The degree of equivalence derived from an RMO key comparison has the same status as 
that derived from a CIPM key comparison.” To further support mutual confidence in the 
validity of calibration and measurement certificates issued by participating institutes so-called 
supplementary comparisons shall be carried out in addition. The metrology community has 
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followed these documented requirements to carry out international comparison measurements 
13 on a regular basis, whose results are continuously registered and made available to the 
public in the BIPM key comparison database (KCDB) 13 . 

The issue of linking by means of joint participants is already addressed in the 
MRA, however without prescribing a particular mathematical linking procedure. Different 
approaches for robust statistical linking of the results of comparisons, e. g. by transferring 
a key comparison reference value (KCRV) from one comparison to the other, have been 
discussed and proposed in literature and have been applied for the analysis of comparison 
results, too 01, ©, Q, El. 

Most of the proposed approaches for robust statistical linking follow the basic 
assumption of the MRA, namely that the results of an RMO key comparison have to be linked 
to the results of a CIPM key comparison. The CIPM key comparison in this concept generally 
is regarded as a primary comparison, both with respect to the chronological order and with 
respect to the participating laboratories, because “... participation in a CIPM key comparison 
is open to laboratories having the highest technical competence and experience, normally the 
member laboratories of the appropriate Consultative Committee.” (see sect. 6.1 of ||T1)- In 
this sense the linking of CIPM key comparisons and RMO key comparisons can be regarded as 
a hierarchical process. However, in 0) approaches for linking of comparisons through global 
minimisation by application of generalised least squares methods El were shortly discussed 
too, which can be regarded as examples of non-hierarchical linking methods. 

In the field of dimensional metrology which is under the responsibility of the 
CIPM Consultative Committee for Length (CCL), international comparisons were started 
immediately after the MRA was signed. The results of the first CIPM key comparison on 
gauge blocks e. g. were already published in 2002 m. 

Based on the experience gained from the international comparisons in the field of 
dimensional metrology, the CCL working group on dimensional metrology (CCL-WGDlvjil) 
has published a document which describes the characteristics of dimensional comparisons 
and which also proposes an alternative comparison scheme, the so-called CCL-RMO key 
comparison scheme Eol. The background and the details of the proposed comparison scheme 
are not repeated here, however some core statements from this document are cited; 

• The CCL-RMO key comparisons follow the idea of several comparisons mutually linked 
together, without the necessity of a top level comparison delivering a KCRV. 

• The classical scheme of CCL and RMO key comparisons can be considered to be a 
special case of the more general CCL-RMO scheme. 

• In terms of linking laboratories within and across regions by calculation of their 
respective degrees of equivalence, the classical (hierarchical) and the CCL-RMO 
comparison scheme can be regarded to be equivalent. The linking of the comparisons 
is guaranteed by common participation of selected laboratories. 

InlfTTl a generalised formalism was presented for linking a CIPM key comparison with 
similar regional key comparisons. As an example, this procedure was applied to link the 
results of a regional (SIM) gauge block comparison llT2ll to the results of the hrst CIPM 
key comparison on gauge blocks ||9l. In the applied procedure, a so-called linking invariant 
parameter is calculated based on the results of the linking laboratories, i. e. those laboratories 
which participated in both comparisons, which in this case used different transfer standards. 

This proposed method was discussed and well accepted in the CCL-WGDM, however it 
was also criticised that it still required to choose one comparison as primary. This condition 

$ In June 2009 the working group structure of CCL has changed: MRA related work is now dealt with by the 
WG-MRA while linking issues are addressed by the task group on on KC linking (TG-L) 
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is not in full alignment of the linking procedure to be applied for the analysis of comparisons 
according to the proposed CCL-RMO key comparison scheme. It was argued, that an adapted 
‘distributed linking’ procedure should be developed which should provide a statistically 
rigorous linking taking into account all the available information of the comparisons to be 
linked, however without the necessity to choose one comparison as a primary one. 

In this paper a Bayesian approach for implementation of such a ‘distributed linking’ is 
presented and an example for linking of two gauge block comparisons is discussed, which 
refers to the data published in j^. lfTSI . 

2. The scenario 

Let us assume the following scenario; A group of laboratories (group A) has participated in 
a key comparison of a travelling standard A and another group of laboratories (group B) has 
participated in a key comparison of a travelling standard B. The task of each group was to 
measure one single measurand of the respective travelling standard (for example the length 
of a gauge block). In order to link the results of the two key comparisons, a certain group 
of laboratories (group C) has participated in both key comparisons, i. e. the laboratories of 
this group have been members of group A as well as of group B and thus have measured 
both travelling standards. All laboratories have reported the respective measurement results 
and the standard uncertainties associated with them. In addition, the laboratories of group C 
have reported the covariances associated with the measurement results of the two different 
travelling standards. It should further be assumed, that the standards have been stable, i. e. 
that they have not changed their properties in terms of time during the key comparisons. 

In order to simplify the organisation of the data, we assume a labelling of the laboratories 
by assigning, independently of the membership in one of the groups, a unique natural number 
to each of them, starting by the number one in ascending order. In the following we will use 
the labels as indices for the respective data. It should be understood, although we are using an 
ordered set of indices, it is totally arbitrary which index is assigned to which laboratory. Any 
permutation of the indices used would also serve. 

After having assigned the indices, we are able to refer to the index sets of the laboratories, 
rather than to the groups. Those laboratories, which have e. g. only measured the travelling 
standard A are within the index set Ta \ 3b- 

3. The information available 

After having assigned labels (indices) to the laboratories, we are able to summarise the 
available data in a clear way just by using the same indices for the laboratories and the data. 
In the following, we will denote the measurement results concerning the travelling standard 
A obtained by the laboratories of group A by xaj {i G 3a) and the standard uncertainties 
associated with this results by u{xA.i) {i G 3a). Accordingly, we will denote the measurement 
results concerning the travelling standard B obtained by the laboratories of group B by xb,; 
{i G 3b) and the standard uncertainties associated with this results by m(xb,j) {i G 3b). It then 
only remains to deal with the covariances, associated with the measurement results concerning 
both travelling standards, as reported by the laboratories of the group C. In the following we 
will denote these covariances by u{xA,i,XBj) {i G (3An3B)). Equivalently we may use the 
correlation coefficients, defined by 


'•ab,; = 


u(xA,i,XB,,j 


i G (3An3B), 


( 1 ) 
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instead of the covariances. We will frequently use the correlation coefficients as abbreviations 
in the subsequent calculations. 

Under the prerequisite, that the two travelling standards under consideration have been 
stable during the key comparison, the assumption is justified, that all measurement results 
are the experimental realisation of only two different quantities, i. e. the measurands of the 
travelling standards A and B, denoted in the following by Fa and Fb, respectively. In addition 
we assume, that the measurement results have already been corrected for possible systematic 
deviations in a suitable way and the uncertainties associated with the corrections have been 
included in the combined standard uncertainties associated with the corrected results. 

Note, that we use upper case letters to denote the measurands and not the respective 
measured value of a measurand, which we will denote by lower case letters throughout the 
text. We will however not differentiate between a measurand and its value (an ideal value to be 
estimated, sometimes called true value) and use the same letter for both, since the respective 
meaning follows from the context. 

4. Establishing key comparison reference values 

4.1. Preliminary remarks 

Using the ordinary arithmetic mean of the measured values to calculate a key comparison 
reference value (KCRV) can not be justified, because in doing so none of the known 
uncertainty values are taken into account and thus we do not make use of all available 
information, as required by the Guide to the Expression of Uncertainty in Measurement 
(GUM) ifTSl . Therefore we shall apply the Bayesian theory of measurement uncertainty lfT4l . 
which is entirely based on Bayesian statistics and the principle of maximum entropy (PME). 
The idea of this theory is first of all to establish a joint probability density function (pdf) of 
the measurands under consideration by taking into account the data obtained by measurement 
or otherwise given, together with their associated uncertainties, as well as physical relations 
and prior information about the measurands. This joint pdf expresses the state of incomplete 
knowledge about the measurands and is subsequently used to calculate the expectations of the 
measurands and its associated covariance matrix. 

4.2. The joint probability density function 

Since the measurements of all laboratories participating in the key comparison are assumed 
to be independent, we can assign probability density functions to the data of each laboratory 
individually. These pdfs belong to two different types only, because the data types of the 
groups A and B, excluding the data of group C, are similar and will thus lead to the same type 
of pdf. The data type of group C, however, is dissimilar from the data types of the groups A 
and B, excluding the data of group C, and hence we have to expect to get a different type of 
pdf for this group. 

In order to assign pdf^to the data, we follow a suggestion of Jaynes ina, Kni, imao 
use the principle of maximum entropy (PME) for this purpose. The PME requires to maximise 
the relative information entropy 1(191 . Il20l by a suitable pdf under the constraints imposed 
on the pdf by the normalisation condition and the data. In the process we make the usual 

§ We will use here and throughout the document the abbreviation 

p{X\Y) = px\y{x\y) 

for the conditional pdf, in order to simplify the notation of the formulae. 
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assumption, that the measured data and the associated variances and covariances, respectively, 
may be equated with the respective expectations. 

The PME yields under the circumstances at hand 


p{xA,i,u{xA,i) I Ya) = 
p{xB,i,u{xB,i) I Tb) = 


1 


u{xA,i)V^ 

1 


exp 


u{xB.i)V^^^'^ \ 2m2(xb_,) 7’ 


{YA-XAjf 

2u^{xA,i) 

{Yb-xbj)^ 


i €3 a\'3b, (2) 
iG3B\'3A, (3) 


and 


p{xA,i,XB,i,uixA,i),u{xB,i),rAB,i I iAjie) = 


271 : m(.xa,!)“ 


X exp — 


2(1-^ab./)L u^ixA,i) 


{YA-XA,i)^ ^ rABAYA-XA,i)iYB-XBj) ^ {YB-XB.if 


u(xA,i)u{xB,i) 


U^{xB,i) \ J 

i S 3a n Jb ■ (4) 


These pdfs can be esteemed as sampling distributions an thus be used to establish the 
likelihood function in the usual way. Note that the values xaj {i G 3a) xb.; {i S 3b) 
as well as their associated uncertainties m(xa,,) and u{xbj), respectively, are given data, while 
Fa and Fb are unknown quantities. 

The prior pdf p(Fa,Fb) would allow us to express our knowledge about the quantities Fa 
and Fb. However, in order to be conservative, we use Jeffreys’ prior US here, i. e. we assign 
a constant to this prior pdf. If more knowledge about Fa and Fb would be available, the prior 
could be changed accordingly. 

Using Bayes’ theorem we obtain the posterior pdf 

/7(Fa,Fb |S))=Ce-^'/^ (5) 

with the set D representing all known data. 


E 


+ E 

/G(3:An3B) 


(Fa-7Ca,;) 

U^{xA,i) 

1 r 


^ V (^B XbjA 
iei3%A) 

(Fa - xajA ^ rAB,i{Y a - xaj) (yn - 

U^{xA,i) u{xA,i)Ax^,i) 


XB,i) 


(^B -XB.iY 

U^{xB.i) 


( 6 ) 


and a normalisation constant C. This is the posterior pdf of the quantities Fa and Fb given the 
data D. 

After some algebraic transformations, equation (|6]l can be written as 


„2 1 ( {Ya-Sa? {YA-yA){YB-yB) iYB-yBf \ 2 

1-^abV «^(fa) m(fa)m(fb) m^(fb) J 


(7) 


with 


FA 


bsi +CS2 
ab — c^ 
csi +as2 
ab — c^ 


( 8 ) 


FB = 


(9) 
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m(ta) = ^ 

/ ^ 

1 ab — c^ 

m(tb) = ^ 

j a 

1 ab — c^ 


c 

''AB = 

s/ab 


( 10 ) 

( 11 ) 

( 12 ) 


2 _ Sa)'^ ^ y-> Ss)^ 

^ ”<e(a7\3B) 

^ y 1 { {xa^-SkY ^ ?'AB,i(^A,;-yA)(-^Bi-yB) ^ (xB,i-yB)7 

iGpfnaB) V «^(aa,,) M(xAi)M(xBi) u^{xQ.i) ) 


and the abbreviations 


Si 


E 

3a \ 

E 


/G(a^aB) “^(-^A.i) ,G(a^ag) (1 '■ab.()“^('^a.() 

1 ^ 1 


ig(3b\3a) 


U^ixBJ 


^ (T 


(l-dB,>^(3^Bi) 


/G(aAnaB)'' AB.i 
'‘AB.i 


c= y _ 

H^AnOj,) (1 - '•ab,X^a.i)m(ab,;) 


XA. 


= (3a\3b) H3%b) ^ " '‘ABi WixA,i) m(xa,')m(xb,') 


aa. 


?‘AB,!AB,( 


(14) 

(15) 

(16) 

(17) 


S2 


E 

z€(3b \^a) 


XB,i 

u'^{xB,i) 


E 

lG(3An3B) 


1 


1 -r^ 

^ '^AB.i 


ab,i 

m^(ab,() 


yAB^i'AA,,- \ 

u{xaj)u{xbj) J 


(18) 


Note, that the auxiliary quantities a, b, c, si, S 2 , and thereby the quantities yA, Sb, u{yA), 
u(yB), Pab as well, depend only on the data obtained by the measurement. Thus does not 
depend on the quantities Fa and Tb and can hence be absorbed by the normalisation constant. 
Therefore, after renormalisation, we finally get 


p(Ta,Tb|®) =- ^= 

27tu(yA)u(yB)^ 1 - rig 

^ I 1 f {YA-yAf (yA-FA)(yB-FB) , (Ib-Fb) 

2(l-r^)v u^iyA) u{yA)u{yB) u^ih) 

i. e. the posterior pdf of the quantities Ta and Yb is a bivariate Gaussian pdf. 



(19) 


4.3. The key comparison reference values 

Since the posterior pdf of the quantities Ta and Tb is a bivariate Gaussian pdf, it follows that 


E [Ta] = yA 


E [Tb] = yB 


( 20 ) 
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Vy 


f fABM(FA)M(FB)\ 

\fABu{yA)uiyB) m ^( fb ) J 


( 21 ) 


their associated covariance matrix. We now identify the values ja and ye as the key 
comparison reference values and the matrix Uy as their associated covariance matrix. 

For convenience we summarise the results, using the covariances rather than the 
correlation coefficients. 


fl= Y, 


1 

U^{xA,i) 


L 

(e(3An3B) 


_ u^{xB,i) _ 

U^ixAjWixBj) - U^{xA,i,XBj) ’ 


( 22 ) 


E 


U^{xA,i 


'6(3b\3a) <e(3l^aB) «^(^A,')«^(^B,-) - lfl{xA,UXB,i) ’ 

m(xa,aab,/) 


Si = 


S2 = 


XA.i 


_ U^{xB.i)xA,i - u{xA,i,XB,i)xB,i 

E(3t\3B) ^ie(a7;^aB) «^(AAi)M2(xB,') -m2(xa,,-,Xb,-) 

^ M^(xa.,)xb,( - u{xA,i,XB.i)XA,i 


Ab,, 


<-e(3t\aA) iG(3l^3B) «^(^A,')m2(xb,,-) - M2(xA,',XBi) 

Key comparison reference values (KCRVs): 

bsi +CS2 


FA = 


Fb = 


ab — c^ ’ 

CSl +flS2 
ab — c 


2 ■ 


Standard uncertainties of the KCRVs: 

m(fa) = 


ab — c^ 


(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 


Covariance of the KCRVs: 


5. Degrees of equivalence 

The degrees of equivalence (DOE) are defined to be 

dA,i = XA.i - FA , i G 2 a , (32) 


and 


dB.i = Ab,( — Fb , 


i G 2 b . 


(33) 
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These values give the deviations of each laboratory from the respective KCRV. 
Following the rules for the propagation of variances, we obtain 


u{dA.i) = 1 

Ju^ixhd)- 

u^iyA) , 

/ £ Ta , 

(34) 

m ( c / b .() = 1 

Ju^{xQ.i) - 


/ £ Tb j 

(35) 


for the respective uncertainties. The negative sign under the square root is due to the 
correlations u{yp^^x^±) = and u{y-g,,x^ k) = as shown in the appendix. 

6. Conformity test 

After the evaluation procedure has been accomplished, we have to check, whether the 
determined results do conform to the data given. This can be achieved in the following way. 
We observe, that on the one hand is represented by equation (|6l) and on the other hand by 
equation (|7]i. Thus, we can calculate the expectation E [x^] in two different ways, yielding 


E[X^] =N = 2 + q\ 

(36) 

with 


N = cardTA + card Ob 

(37) 

and q^ given by equation (fTsTl. Thus, we obtain 


q^=N-2. 

(38) 


Note that the right hand side of this equation may be regarded as the number of the degrees 
of freedom, because it is just the total number of measured data diminished by the number of 
estimated quantities. 

The condition ( l38l l can, however, not be expected to be strictly fulfilled, but A — 2 is the 
most probable value of q^, because if the data are actually statistically distributed according 
to the model and the measurement uncertainties given, equation dMt is the consequence. A 
significant deviation from this result is a signal that either the model is wrong or the data are 
suspect. Assuming the model to be correct, an outcome > N — 2 leads to the conclusion, 
that there is a strong possibility that either some of the measurement uncertainties have 
been underestimated or the respective measurement results contain uncorrected (or unknown) 
systematic deviations. Thus 

q^<N-2 (39) 

can be taken as an indication of the conformity of the estimated results with the data. 

7. Examples 

In order to demonstrate the application of the method proposed in this paper, we show two 
examples. Our first example is a synthetic one. Data have been produced by a simulated 
sampling from Gaussian distributions as given in equations ^ to (|lli, using the parameter^ 
Yp^ — \ 10, Oa = 20, Tb = 120, Ob = 50, p = 0,5 and a sample size of n = 50. Subsequently 
the sampled data have been evaluated by using the usual formulae for the sample mean, the 
sample variance and the sample covariance. The resulting data thus obtained as well as their 
evaluation are shown in table [T] and table |2] 




A Bayesian approach to the Unking of key comparisons 


9 


Table 1. Data of the synthetic example obtained from the sampled data by application of the 
usual formulae for the sample mean, the sample variance and the sample covariance. 


Institute 

XA.i 

u{xA,i) 

XB.i 

u{xBj) 

^AB,i 

LAB-01 

113,4 

2,9 




LAB-02 

112,1 

2,8 




LAB-03 

113,0 

2,5 




LAB-04 

110,6 

2,6 




LAB-05 

109,4 

2,4 




LAB-06 

107,0 

2,6 




LAB-07 

104,7 

2,8 




LAB-08 

109,0 

2,6 




LAB-09 

111,0 

2,4 

120,1 

6,5 

0,8 

LAB-10 

109,4 

2,8 

117,3 

7,3 

0,8 

LAB-11 

111,1 

2,8 

125,0 

6,4 

0,8 

LAB-12 

115,3 

2,4 

135,7 

6,7 

0,7 

LAB-13 



129,7 

6,1 


LAB-14 



129,1 

7,5 


LAB-15 



125,0 

7,1 


LAB-16 



123,6 

6,6 


LAB-17 



123,0 

6,9 



As can be seen, the data do pass the conformity test. However, the KCRVs are slightly 
larger than the values used for the simulation. This is caused by the fact, that the correlation 
coefficients estimated from the sampled data tend to be too large compared to the correlation 
coefficient used for the simulation. This demonstrates the effect of small sample sizes. 

As a second example data from the CCL-Kl il and SIM.L-Kl fUl gauge block 
comparisons have been evaluated. The three institutes CENAM, NIST and NRC participated 
in both comparisons and hence serve as linking laboratories. For convenience, the results 
and their associated standard uncertainties, as reported in the respective publications, are 
repeated in table[3 Values for the covariances have unfortunately not been reported, although 
systematic deviations present when measuring the gauge blocks in each of the linking 
laboratories with the same measuring system inevitably cause correlations, even if a correction 
for these systematic effects, as usually good practise, has been ^plied. Thus, following the 
recommendation given in the GUM in lack of this informatioifll all covariances have been 
assumed to be zero, knowing that this is certainly not true. 

The degrees of equivalence for the deviation from the nominal values and their associated 
standard uncertainties, as well as the key comparison reference values and their associated 
standard uncertainties for the groups A and B, respectively, obtained by the application of the 
linking method proposed in this paper are given in table |4] for the original data. It emerges, 
however, that the data do not pass the conformity test. 

Assuming that the evaluation model is correct and that the respective measurement 
results have been corrected for systematic deviations, the test result leads to the conclusion, 
that there is a strong possibility that some of the measurement uncertainties have been 

II a and p denote the standard deviation and the correlation coefficient, respectively, of a Gaussian distribution. 

It should be noted, that this recommendation is consistent with the modus operandi of the Bayesian theory, where 
only known information is taken into account. 
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Table 2. Degrees of equivalence for the deviation from the nominal values and their associated 
standard uncertainties, as well as the key comparison reference values and their associated 
standard uncertainties for the groups A and B, respectively, obtained by the application of the 
linking method proposed in this paper for the synthetic example. 


Institute 

dA.i 

u{dA.i) 

dsj 

u{dBj) 

LAB-01 

2,491 

2,815 



LAB-02 

1,191 

2,712 



LAB-03 

2,091 

2,401 



LAB-04 

-0,309 

2,505 



LAB-05 

-1,509 

2,296 



LAB-06 

-3,909 

2,505 



LAB-07 

-6,209 

2,712 



LAB-08 

-1,909 

2,505 



LAB-09 

0,091 

2,296 

-3,779 

6,196 

LAB-10 

-1,509 

2,712 

-6,579 

7,030 

LAB-11 

0,191 

2,712 

1,121 

6,091 

LAB-12 

4,391 

2,296 

11,821 

6,405 

LAB-13 



5,821 

5,775 

LAB-14 



5,221 

7,238 

LAB-15 



1,121 

6,822 

LAB-16 



-0,279 

6,300 

LAB-17 



-0,879 

6,614 


Ta = 110,909, uIja) = 0,698 


ye = 123,879, M(yB) = 1,966 


q^iN- 

2) =0,89 



underestimated. Inspecting the uncertainty values given in table |3 it is obvious that 
.. INMETROl submitted very optimistic uncertainty claims relative to conventional 
capabilities” ifT^ (the corresponding value is marked in bold face in the table). Since the 
gauge blocks have been measured at INMETRO by two completely different instruments, 
involving different staff members (data denoted by INMETROl represent results from a 
research grade instrument, while those denoted by INMETR02 represent results from the 
routine gauge block calibration service offered by INMETRO) HU, it is reasonable to 
question the very low uncertainty stated for INMETROl. Therefore, this uncertainty value 
has been increased here from 4,0 nm to 11,2 nm, which is the smallest value in order that the 
data pass the conformity test. The respective results obtained with this change are given in 
table |5] 

A comparison of the results given in table |4] with those given in table |5] shows, that the 
reference value for group B and its associated standard uncertainty has noticeably changed. 
This strongly emphasises, that it is imperative to do a rigorous conformity test. The E„ 
criterion, for example, did not reveal any problem with the data from INMETROl during 
the SIM.l-Kl comparison. In fact the values \Eg$ \ as well as jEj for the evaluation of the 100 
mm steel gauge block data have even been comparable to those for the NRC lfT2l . The usual 
x'^ criterion at a 5 % confidence level also fails to indicate a problem with the data. 
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Table 3. Results for the deviations from the nominal length and their associated standai'd 
uncertainties as taken from the CCL-Kl 0 and SIM.L-Kl (TU gauge block comparison 
reports (steel gauge block, nominal length 100 mm). 


Institute 

■*A.,/nm 

“(■*A,r)/nm 

jCB,;/nm 


METAS 

-96,0 

13,0 



NPL 

-140,0 

33,0 



BNM-LNE 

-110,0 

16,0 



KRISS 

-104,3 

20,6 



NRLM 

-89,4 

16,3 



VNIIM 

-104,0 

15,0 



CSIRO 

-114,0 

16,0 



NIM 

-90,0 

10,3 



NIST 

-117,0 

17,9 

-100,0 

18,0 

CENAM 

-119,0 

18,7 

-93,0 

23,0 

NRC 

-126,0 

24,0 

-124,0 

26,0 

INMETRO 1 



-98,0 

4,0 

INMETR02 



-68,0 

29,0 

INTI 



-104,0 

21,0 

CEM 



-148,0 

17,0 


Table 4. Degrees of equivalence for the deviation from the nominal values and their associated 
standard uncertainties, as well as the key comparison reference values and their associated 
standard uncertainties for the groups A and B, respectively, obtained by the application of the 
linking method proposed in this paper (results obtained for the original data in table[3). 


Institute 

£/A,i/nm u{dA.i)/nm 

iiB,!/nm 

u{dB,i)/nm 

METAS 

7,6 

12,1 



NPL 

-36,4 

32,6 



BNM-LNE 

-6,4 

15,2 



KRISS 

-0,7 

20,0 



NRLM 

14,2 

15,6 



VNIIM 

-0,4 

14,2 



CSIRO 

-10,4 

15,2 



NIM 

13,6 

9,1 



NIST 

-13,4 

17,2 

0,5 

17,6 

CENAM 

-15,4 

18,1 

7,5 

22,7 

NRC 

-22,4 

23,5 

-23,5 

25,7 

INMETRO 1 



2,5 

1,7 

INMETR02 



32,5 

28,8 

INTI 



-3,5 

20,7 

CEM 



-47,5 

16,6 


Sa = -103,6, 

u{yA) = 

4,9 nm 



yB = -100,5, 

uIjb) = 3,6 nm 



qy{N-2) = 

., 07 (conformity test failed) 
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Table 5. Degrees of equivalence for the deviation from the nominal values and their associated 
standard uncertainties, as well as the key comparison reference values and their associated 
standard uncertainties for the groups A and B, respectively, obtained by the application of the 
linking method proposed in this paper (results obtained after the uncertainty of the institute 
INMETROl has been increased from 4,0 nm to 11,2 nm). 


Institute 

d\j/nm u{dA,i)/nm 

c/B.j/nm u{dBj)/nm 

METAS 


7,6 

12,1 



NPL 


-36,4 

32,6 



BNM-LNE 


-6,4 

15,2 



KRISS 


-0,7 

20,0 



NRLM 


14,2 

15,6 



VNIIM 


-0,4 

14,2 



CSIRO 


-10,4 

15,2 



NIM 


13,6 

9,1 



NIST 


-13,4 

17,2 

6,7 

16,6 

CENAM 


-15,4 

18,1 

13,7 

22,0 

NRC 


-22,4 

23,5 

-17,3 

25,1 

INMETROl 




8,7 

8,9 

INMETR02 




38,7 

28,2 

INTI 




2,7 

19,9 

CEM 




-41,3 

15,6 

Sa = -103,6, M(yA) =4,9 nm 
ye = —106,7, u{yB) =6,8 nm 
q^/{N-2) = 1,00 


8. Conclusion 

In this contribution we have derived formulae for the linking parameters key comparison 
reference value (KCRV) and degree of equivalence (DOE) for analysis of two comparisons 
based on Bayesian statistics, i. e. taking all known information into account but without 
the necessity to choose one comparison as a primary one. The applicability of the linking 
approach was demonstrated on an example from gauge block metrology. 

The proposed approach can be described as a type of ‘distributed linking’, which has 
been looked for especially in the area of dimensional metrology within the working groups of 
the CCL. 

This ‘distributed linking’ method can, for example, be used to compare the results of two 
comparisons which were started in parallel using different transfer standards and where it does 
not seem meaningful to choose one comparison to be primary (defining a KCRV) and the other 
secondary (linking its comparison results to the KCRV of the first). Such comparisons with 
more than one loop have been organised if the number of interested participants became rather 
large, like e. g. in the finished line scale comparison EUROMET.L-K7 with 31 participants 

EH. 

In a strict Bayesian sense, using the ‘distributed linking’ approach takes all known 
information from both comparisons into account to determine the linking parameters, i. e. 
all the measurement results and their uncertainties. Instead, in application of the hierarchical 
linking approaches one could in principle argue that not all of the information is taken into 
account, because e. g. the results gained from a RMO key comparison are not taken into 
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account to reflect or re-analyse the results of a CIPM key comparison, although by the 
twofold involvement of the linking laboratories the level of knowledge is increased for both 
comparisons. On the other hand the results of a subsequent RMO comparison linked to a key 
comparison by the described ‘distributed linking’ approach would also change the KCRV of 
the (already finished) key comparison. We see that this would impose new challenges for the 
operation of the KCDB. 

The described approach can be extended to link more than two comparisons. Moreover, 
additional knowledge, as e. g. results from prior measurements, may be taken into account by 
choosing a suitable prior pdf. 


9. Acknowledgement 

The authors gratefully acknowledge the valuable comments of A. Balsamo from INRIM 
(Italy) and the suggestion to add the information given in appendix 2. 


Appendix 1 (Calculation of covariances) 


In equations (l34l i and dTSl l the relations u{yj^,XQk) — and u{ys,XBii) = u^iys), 

respectively, have implicitly been used. In order to show that they are valid, we introduce 


the estimators 


bS 1+ cS2 
ab — c^ 


(40) 


and 


Yb 


cS 1 aS 2 
ab — c^ ’ 


(41) 


where 


and 


5i 


L 

'6(5a\3b) 


U^{xA,i) 


S2= E 

ig(3b\3'a) 


m2(xb,,) 


!G(a7n3B) «^(-^A,,')M^(jCB,i) -«^(-^A,/,-«B,,-) 

U^{xA,i)XB,i - u{xA,i,XB.i)XA,i 
iepfnJg) u^(xa.,)u^(xbj) - u^(xA,i,XB,/) ' 


(42) 


(43) 


are functions of the random variables Aa,/ (i € 3a) and Ab,, (i € 3b)- We assume the measured 
values XA.;- (i G 3a) and xb./ (i G 3b) to be a realisation of these random variables. 

Since E[Aa,,] = Ya, E[Ab,,] = Tb, Var[AA,i] = u^(xA,i), Var[AB,,] = u^ixBj) and 
Cov [Aa,,-, Ab,,] = m(xa,,',xb,,) (i S 3a, i G 3b) it can be verified, that 


E [5i] =fl}A-ciB, E [52] = Wb - cTa , (44) 


as well as 

Var[ 5 i]=a, Var[ 52 ]=f 7 , Cov [51,^2] =—c ( 45 ) 

is valid. Using these results, we obtain 

E [Ta] =yA, E [Tb] = JB, (46) 

Var[yA] =M(yA), Var[TB] =M(i'B), Cov [Taj^b] = M(yA,yB), (47) 

i. e. the estimators Ta and Tb yield the correct results as given by equations (IZTT i to dul l. 
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COV [5l,XA,i:] = COV = 1 

(48) 

and 

CoV [5l,7fB,(t] = CoV [52,7fA,t:] =0. 

(49) 

Thus 

u{yA,XA,k) = Cov [TA,XA,i:] = u^{yA) 

(50) 

and 

u{yB,XB.k) = Cov [TB,7fB./t] = k/{yB) ■ 

(51) 

This is what we wanted to show. 


Appendix 2 (Uncorrelated measurement results) 


Here we consider the possibility that the measurement results of the linking laboratories 
(group C) have to be treated as if they were uncorrelated, because the respective covariances 
have not been reported. In this case we have to set all covariances M(xA,,,XB,i) (i G (TaD^b)) 
equal to zero in the equations (|22]) to (l26l). This yields 


(52) 


(53) 

C = 0, 

(54) 

^A,i 

(55) 

4 ’ 

(56) 


Thus, we obtain from equations (l27l) and (l28l) for the key comparison reference values 

S2 


■Si 

yA = — 
a 


and 




and from equations (|2^ and (l30l l for their associated standard uncertainties 


i{yh) = 


s/a 


and 


<(yB) = 


s/b' 


(57) 


(58) 


i. e. the usual weighted mean results of for the KCRVs separately for each group participating 
in the key comparisons. The results of the groups are independent of each other as indicated 
by the covariance of the KCRVs, u{yf^,y-a) = 0, as obtained from equation OTI) . This 
demonstrates that the linking procedure essentially depends on the covariances reported by 
the linking laboratories. 
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